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ABSTRACT 

Detailed heat flow studies were carried out with 
bottom-hole temperature data, geological formation data and 
measured thermal conductivity values obtained from core 
Samples taken from wells in west central Alberta. 

The results indicate that high vertical geothermal 
Qradients are observed bordering the eastern limit of the 
disturbed belt. Also, detailed subsurface temperature 
analysis indicates that lateral temperature gradients exist. 
These gradients originate along the disturbed belt and 
strike northward and northeastward. They are probably 
associated with Pluid motion within the permeable 
Cretaceous, Mississippian and upper Devonian formations. 

Calculated conductive heat flow was found to increase 
with depth. Heat flow across the Mesozoic formations was 
found to be 69.19mWm~?, while 95.03mWm~? was obtained within 
the Paleozoic sediments. The difference in heat flow values 
above and below the Paleozoic erosional surface is 
attributed to fluid motion. A hydrodynamic model can explain 
the thermal regime. The model suggests that gravity induced 
water flow in the Rocky Mountains and the foothills region 
and along permeable Cretaceous, Mississippian and Devonian 
Strata is a major influence on the thermal regime. 

Finally, the evidence suggests that steady-state 
"background" heat flow can only be determined from careful 
measurements below the sedimentary strata in the deeply 


buried Precambrian basement formations which may be less 
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1. INTRODUCTION 


1.1 GENERAL 

Heat flow from the deep interior of the Earth toward 
its surface provides the energy to drive tectonic processes 
and to cause metamorphic and igneous activity. A study of 
snchemheatartlowiwand related thermaleprocesses wisi therefore 
important to better understand the nature of the Earth. 

For a long time it has been known that temperature 
increases with depth into the Earth, and this temperature 
gradient produces an outward heat flow. The main source of 
the heat which flows outward from the interior of the Earth 
is believed to be from radioactive decay of long-lived 
isotopes (Birch et al, 1968; Roy et al, 1968) supplemented 
by slow cooling of the Earth as a whole (Bott, 1982). 

Although a small fraction of the heat flowing outward 
is transmitted to the Earth's surface and lost to the 
Surroundings, the majority 1S converted to other forms of 
energy to drive tectonic processes and to cauSe igneous and 
metamorphic activity. The near-surface heat flow ranges from 
0-200 mWm-2 over 99 percent of the Earth's surface (Elder, 
1964), land+aBrecentocompilationi by SSelaters etyi adrce (1980) 
indicates an average worldwide value of 82mWm~’. In 
addition, there are localized thermal areas where the heat 
fLUxMenSArUp Btommd0 timesfthemnormad vyaluey, andgthese areas 


are characterized by hot springs and steaming ground. 
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Heat may be transmitted to the Earth's surface by 
conduction through the solid crystalline rocks, by water 
motion; sand by’ upward flowing volcanic material: Heat 
conduction occurs throughout the whole Earth, while volcanic 
action occurs in limited areas and may provide information 
from depths as great as 100 km (Elder, 1981). Hydrothermal 
systems are much more widespread, but can provide 
information only to depths of about 10 km. 

Thermally conducted heat flow per unit area is equal to 
the temperature gradient (VT) multiplied by the thermal 
conductivity K: Q=VTxK. Applying this to the Earth, the heat 
which escapes by conduction from the interior at a location 
on the surface can be determined by measuring 
(i) the temperature gradient just below the Earth's surface, 
“and 
(ii) the thermal conductivity of the rocks. 

The redistribution of heat by active groundwater flow 
has been the subject of heat flow studies for as long as 
measurements have been made (Van Ostrand, 1934; Bullard, 
1939). It is clear that convecting systems play an important 
role in heat transport in thermal areas (Lachenbruch and 
Sass, 1977; Bodvarson, 1969). Also there iS evidence that 
heat transport by both local and large scale water motion 
strongly affects the thermal regimes of wider areas such as 
sedimentary basins (Majorowicz and Jessop, 1981a,b; Smith 
and Chapman, 1982). AS a consequence, study of the 


temperature regime and related heat flow in a basin can 
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provide important information about the nature of fluid 
motion and thermal interactions there. 

Understandinig the temperature regime and heat flow in 
basins 1S important for evaluating the geothermal energy 
potential (Rybach and Muffer, 1981; Jessop, 1976; BRGM 
report, 1976). Also, although temperature information has 
not been used directly in hydrocarbon exploration, it is 
speculated that oil and gas deposits may be associated with 
temperature anomalies (Roberts, 1979). Water movement 
transports hydrocarbons and also influences the thermal 
regime, so some relationship between hydrocarbon reservoirs 


and temperature distribution can be expected. 
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1.1.1 Review of Basin Heat Flow Studies from Bottom-Hole 
Temperature Values 

Sedimentary basin heat flow studies can assist in 
locating potential regions for low-grade geothermal energy 
recovery as well as provide information related to the 
generation and migration of hydrocarbons within the basin. 
Also, temperature surveying iS becoming important Lor 
locating groundwater (O' Brien, 1970). 

The temperature distribution within the sediments and 
their thermal conductivities need to be known to describe 
the thermal regime of the basin. Also, the porosities and 
permeabilities of the sediments are important parameters in 
assesSing the presence of fluid and the nature of its 
motion. 

Unfortunately, accurate temperature observations are 
seldom available, except perhaps at widely separated 
isolated locations. However, in most basins, petroleum 
exploration has taken place, and Zone temperature 
measurements have been associated with these efforts. 

When service company logging iS carried out in an 
exploration or development well, bottom-hole temperatures 
are uSually taken during each logging run by a maximum 
reading thermometer as a matter of course. These temperature 
measurements abound with both human and mechanical errors, 
and conclusions from them must be considered with these 
errors in mind. Also, well temperatures are normally 
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and it may take a long time for temperature equilibrium to 
return to the well and its surroundings. This effect must 
also be taken into account when such measurements are used 
for subsurface temperature estimates. 

Although many difficuities and uncertainties occur when 
bottom-hole temperature (BHT) values are used, it must be 
recognized that these measurements are the only data that 
exist in many areas. If many BHT values are used, and 
interpretation based on them made with care, useful 
information should be gained. 

Bottom-hole temperature values from well logs have been 
used previously for basin studies in Europe, North America 
and elsewhere. Harper (1971) and Evans and Coleman (1979) 
were the first to use BHT data in estimating heat flow for 
the North Sea sedimentary basin. Recently Oxburgh and 
Andrews-Speed (1981) used BHT together with estimated rock 
thermal conductivities from meaSurements on 230 samples to 
evaluate heat flow in the south western North Sea. They 
found average heat flow values of 70-75mWm~? which were 
generally higher than the average value of 60mWm~? found on 
the mainland United Kingdom. Also, the heat flow 
distribution pattern exhibited a correlation with the main 
geological features. 

Balling et al, (1980) used BHT to investigate the 
geothermal energy potential of the Danish sub-basin which is 
part of the west European sedimentary basin, and Gable 


(1979) used BHT data to eStimate heat flow in the Paris 
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basin in France. 

Reiter et al, (1982) have also used BHT to study the 
Northern Chihuahua geothermal field in Mexico. They used 
many uncorrected BHT values and have been able to delineate 
areas with heat flows ranging from 60mWm-? to approximately 
90mWm-* in the region. 

Reiter and Jessop (1983) have used BHT to identify heat 
flow regimes in the offshore eastern Canadian basin, and 
Anglin and Beck (1965) and Majorowicz and Jessop (1981a) 
have made regional heat flow estimates from BHT and 
temperatures from shut-in wells in the central Canadian 
plains region. In both the eastern Canadian offshore work of 
Reiter and Jessop (1983) and the plains work of Majorowicz 
and Jessop (1981a), it was suggested that fluid migration 
along permeable strata or along faults and fractures 
strongly affects the heat flow regime. Recently, Yorath and 
Hyndman (1983) estimated heat flow for the Queen Charlotte 
Basin in Canada using BHT data and effective thermal 
conductivity estimates based on the relationship between 
SOmic velocity and thermal conductivity of rocks. Carvalho 
et al, (1980) estimated heat flow using BHT from Sumatra 
basin wells. 

A variety of factors cause variable heat flow patterns 
in sedimentary basins. The basin age, distribution of 
radioactive heat sources, heat conductivity variations, 
magmatic heat sources and thermal effects of groundwater 


flows are among the most important factors influencing the 
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1.1.2 Previous Heat Flow Studies in West Central Alberta 

The first heat flow study in the western Canadian 
sedimentary basin in Alberta was undertaken by Garland and 
Lennox (1962). They measured equilibrium temperatures in two 
wells in Redwater and Leduc and used heat conductivity 
values from samples from nearby wells to determine heat flow 
values of 61mWm-? and 67mWm~? which are rather higher’ than 
the average value for sedimentary basins with Precambrian 
basements. 

Anglin and Beck (1965) were the first to use 
temperature data collected from orl companies. From 
geothermal gradient estimates for 18 sites and assuming that 
the conductivity 1s reasonably constant for the basin, they 
found a northerly trending increase in heat flow values 
which they suggested may be due to crustal thickening. Also, 
they suggested a correlation of enhanced temperature with 
granitic basement intrusions and, furthermore, their results 
showed temperature gradient values that increased eastward 
in the direction of decreasing sediment thickness. 

As mentioned in the previous section, Majorowicz and 
Jessop (1981a) made estimates of regional geothermal 
gradients and heat flow from available temperature data 
(AAPG, 1976) for the interior plains. They found that the 
heat flow values are generally higher than those found in 
other Precambrian platform areas. They found that 
radioactive heat generation estimated on the basis of 


uranium, thorium, and potassium studies by Burwash and 
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Cumming (1976) and Burwash (1979) could not explain the heat 
flow pattern and suggested that fluid motion is a major 
contributing factor to the regional heat flow pattern in the 
region, 

Lam et al, (1982) used BHT values from west central 
Alberta and found a high geothermal gradient anomaly in the 
Hinton/Edson area. They suggested that this anomaly may be 
the result of hot water motion northeastward and upward 
toward the surface from the disturbed region to the west 
along southwestward dipping faults in the Devonian and 
Cambrian sediments. 

A closely spaced magnetic variometer array study 
conducted by Ingham et al, (1983) has indicated a deep 
electrical conductivity anomaly coincident with the 
geothermal anomaly ieyasimaoe by Lam et al, (1982). Ingham et 
al, (1983) suggested an alternate explanation for the 
Hinton/Edson geothermal anomaly in that it may be due to a 
region of partial melt with its top some 5-10 km below the 
surface and extending to considerable depth. However, 
further analysis of BHT data by Majorowicz et al, (1983a) 
and Majorowicz et al, (1983b) Shows that the location “and 
general shape of the geothermal anomaly varies with depth, 
which is inconsistent with an anomaly associated with a deep 
source of partial melt. 

The present work is a detailed study of BHT data from 
the Hinton/Edson area in order to gain better knowledge of 


the temperature regime there. Also, detailed heat flow 
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estimates are made in this work. 


1.1.3 Geology of the Area 

The area of interest lies between latitude 52.7 degrees 
north to latitude 54.0 degrees north and between longitude 
115.8 degrees west to longitude 118.0 degrees west, and is 
part of the interior plains of western Canada. 

Western Canada may be divided into three geological and 
physiological provinces. From west to east these are 
(i) the Cordillera 
Ga jmesthe intervorsplains 
(iii) the Canadian shield. 

The Cordillera has extreme topographic relief and forms 
the western boundary of the interior plains. The rocks are 
mostly miogeoclinal and were deformed during the Cretaceous 
and early Cenozoic orogenies. 

The interior plains, or western Canadian sedimentary 
basin etsmaa = broad @gently Tolling plaine frome the: Rocky 
Mountain foothills in the west to the Precambrian shield of 
Pentel Canada. The thickness of the sedimentary column is 
about 7,000m in the west and gradually thins until it 
di@appears some J000km or ~So ©to the east. The rocks are 
Paleozoic, Mesozoic, and Tertiary. The Mesozoic 1S separated 
unconformably from the Paleozoic. The middle Cambrian 
extends far out from the mountains across the basin. The 
Devonian rocks are reasonably continuous’ throughout the 
basin, and are terminated in the easterly and northeasterly 
directions by post-Paleozoic erosion. The Paleozoic finally 
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bordering the shield. 

The Cretaceous rocks cover the entire basin, and dip in 
south west direction. The Paleozoic, Triassic and Jurassic 
strata are basically marine. The Cretaceous is both 
continental and marine, while the Cenozoic is completely 
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2. DATA COLLECTION AND ANALYSIS 


2.1 Sources of Data 

Three sources of data were used for this work: 
(i) Bottom-hole temperatures (BHT ) from petroleum 
exploration well logs kept at the Alberta Energy Resources 
Conservation Board in Calgary. These temperatures are 
recorded at well bottoms by maximum recording thermometers 
attached to logging equipment. These data are taken from 
about eight hundred wells whose locations’ range from 
latitude 52.7 degrees north to latitude 54.0 degrees north 
and from longitude 115.8 degrees west to longitude 118.0 
degrees west. On the township-range grid, this area includes 
townships 4350 57 and@ranges 13 to 27, west] of the ‘Sth 
meridian. The entire area 1S about 21,072 square kilometers, 
and is shown on the map of Fig. 2.1 where the plotted 
positions indicate well locations. The temperature data have 
not been corrected for effects associated with fluid 
Gereulats on durang drilling. 
(2) The second set of data is taken from the Alberta Energy 
Resources Conservation Board tops file. It consists of lists 
Of geological formation tops with theiriswespective depths 
for each well. The ages of these formations range from 
Tertiary to Precambrian period, with about 95% occurring in 
Mesozoic to Paleozoic strata. The wells used were the same 
as those from which BHT dataswere taken. 


(3) The third set of data consists of thermal conductivity 
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Figure 2.1 The plotted positions on this map indicate the 


locations of the wells from which BHT data were taken. 
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measurements made from core samples. The samples were’ taken 
from industry acquired cores from wells in south central 
Alberta, and the measurements were performed on a _ divided 
bar apparatus by the heat flow group at the University of 
Plbertar, In addi e1onecommthe: sthermalwuconductivitys values, 
anformation on «density, Mipowosity and rock type hasbeen 
Recorded | N1Gunem ee ZesnOws ehes posit yvonseoue the) well sa gtrom 


which these data came. 
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Figure 2.2 The plotted positions show the locations of wells 


were taken DO conductivity 


from which core samples 


measurements. 
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2.1.1 Mapping of Geological Surfaces 

Analysis of the three data sets began by mapping 
important geological surfaces. Three surfaces were selected: 
(Ayeetopeot thes Colorado formation, 

(b) the unconformity separating Cretaceous formations 
from the Jurassic-Triassic formations, 
(c) the top of the upper Devonian formation. 

These surfaces were selected because a large number of 
the bottom-hole temperatures occurs in the Cretaceous and 
Devonian formations. The unconformity surface was considered 
because it represents a surface where different lithologies 
are in contact, and “heat ftlow “estimates’? across such, a 
Surface Sshowldyyield instght into whe pattern and naturegor 
the general heat flow of the area. 

To map the surfaces using the geological formation 
depth data, the area was divided into 225 1x1 township/range 
areas. To obtain the representative depths to the top of the 
Colorado formation for each 1x1 area, the depths to this 
formation for all the wells within each area were summed and 
averaged. These averages were assigned to the centres of the 
1x1 township-range areas and from these a contour map was 
drawi@and 1S Shown in Fig." 20" 

The depths to the unconformity surface in each of the 
1x1 township-range areas were determined in a different way 
from those to the top of the Colorado formation. The 
unconformity surface separates formations of the lower 


Cretaceous from those of the JuraSsic-Triassic. or the 
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Figure 2.3 Depths from the surface of the Hacth to, thes Lop 
of Colorado formation (contour values in hundreds of 


meters). The limit of the disturbed belt TSeincicavec. 
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Figure 2.4 Depths (in hundreds of meters) from the surface 
of the Earth to the unconformity surface which separates the 
lower Cretaceous formation from the rest of the Mesozoic 


formations. 
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MissisSippian. Therefore, in each of the 1x1 township/range 
areas, thevdepths tomthe top of the latest formation#oty ethe 
Jurassic were averaged. When no Jurassic formation occurred, 
the most recent formation of the Triassic or Mississippian 
waS considered. The resulting unconformity surface map is 
showmein Fig. 254% 

The third surface, the top of the upper Devonian 
formation, waS mapped in the same way as the top of the 
Colorado formation. In this case the surface was taken to be 
the contact between the Exshaw or the Exshaw Wabamun and the 
Devonian formations. Figure 2.5 illustrates the depth to the 
contact between the Wabamun and the upper Devonian 
formations. 

Observations from these three maps indicate that the 
sediments dip in approximately the southwest direction. 
However, there are not sufficent data in the extreme 
Southwest corner to completely describe how the sediments 
dip as the Rocky Mountain foothills are approached. 

The slopes of the surfaces in the southwest direction 
as far as data exist are: 

(1) the top of the Colorado formation ~12m/km. 
(2) the unconformity surface ~ 16.1m/km. 
(3) the top of the upper Devonian surface = 17.3m/km. 
The contour interval for the three maps is 200m and the 


contour values are in hundreds of meters. 
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Figure 2.5 Depths from the surface of the Earth to the 
contact between the Wabamun_= and the upper Devonian 
formations. (Contour values are given in hundreds of 


meters). 
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2.1.2 Calculation of Geothermal Gradients 
Intorderstoscalculateythe conductives heat flow ein ra 
given direction across a column of rock material, the rate 
of change of temperature with depth (geothermal gradient) is 
required as well as the average thermal conductivity for the 
column. For a homogeneous isotropic material with steady 
heat flow, two temperature meaSurements at different 
positions in the column are needed to define the geothermal 
gradient. The heat flow can then be estimated by taking the 
Product ofetherthermal gradient vandithe rock conductivity: 
ee ae important to realize that a considerable amount 
of scatter’ is introduced into the BHT values due to 
mechanicalf errous Sin Phecordingsnandinalsoi duenston/ftluid 
circulation during drilling. In order to estimate geothermal 
gradients using these data, statistical methods are 
employed. A computer program based on a linear least squares 
fit technique was used to calculate the geothermal gradients 
Erom=temperaturewdepthuplots ofethe:tdata.h»The ection was 
divided into 3x3 township/range areas. Each of these areas 
was moved every two ranges to the west starting from 
townships 43,44,45, ranges 13,14,15 to townships 43,44,45, 
ranges 25,26,27 and then every two townships to the north. 
This moving window procedure was adopted to enhance the 
geothermal gradients and also to average out any abnormal 
temperatures which might otherwise lead to spurious values 
fon-gthe Pealculatedvagradientsse formaany pauhicular SxS 


township/range. 
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By this procedure, average geothermal gradients were 
Calculated for each 3x3 township/range area for the entire 
area. Figures 2.6-2.8 show temperature/depth plots for three 
3x3 township/range areas. The data represented in Fig. 2.6 
are from an area which falls within the Hinton/Edson 
geothermal gradient anomaly reported by Lam et al, (1982). 
Figure 2.7 is taken from the area covering townships 45,46 
and 47.and ranges 19,20 and 21 and falls within the 
disturbed belt. Figure 2.8 gives data from townships 55,56 
andor) and ranges 15,516 and 1% which “are located® within 
the northeast region of the study area. The geothermal 
gradients obtained in this manner were assigned to _ the 
Centres) Of seach fof the corresponding 3x3 township/nange 
areas for contouring. Figure 2.9 is the resulting geothermal 
Geadient = contour “map. The contour interval sfor this=map is 
Ze C/ kine 

Even though the Earth's heat source is\.fairly constant 
and the heat flow is upward, the vertical heat flow varies 
both laterally and vertically. In sedimentary basins, fluid 
motions and rock inhomogeneities are the main causes’ for 
Such variations. It is therefore necessary to calculate 
geothermal gradients over various depth intervals, the same 
way as is done aye sonic velocities in seismic 
interpretation. Four lithological sections are considered 
here. Their surfaces correspond to those mapped in Figs. 
2.3-2.5. These intervals include formations from 


(1) the top of the Colorado formation to the unconformity 
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Figure 2.6 Temperature/depth plot for townships 581,52,53 and 
ranges 19,20,21 showing the lineeter Leastes SQUarTeS. sutited 


line. The broken lines define the 95% confidence limits. 
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Figure 2.7 Temperature/depth plot for townships 45,46,47 and 
Parcecmmlgn70r 2 ue SnoOwWlNG set ncuml ied meleast SQucresatatted 


line. The broken lines define the 95% confidence limits. 
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Figure 2.8 Temperature/depth plot fOnMCOWN Sha Domo, O0,07 sand 
ranges 15,16,17 showing the linear least squares fitted 


line. The broken lines define the 95% confidence limits. 
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Figure 2.9 Geothermal gradient contour map from the 3x3 


township/range moving window estimates. 
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surface, 

(2) the top of the Colorado formation to the upper Devonian 

Surface, 

(3) the surface of the unconformity to the Precambrian 
basement, 
(4) the upper Devonian surface to the Precambrian basement. 

Using the same linear least squares technique as 
before, geothermal gradients were calculated for each of 
these intervals. In all cases the 'spread' was calculated as 
well. It is the standard error of estimate associated with 
the least squares fitted line. 

Figures 2.10 and 2.11 show temperature/depth plots and 
calculated geothermal gradients for the 3x3 areas given by 
BOWNSRIDS Sl, ozs and 53° “and: “ranges” 197% 20 and Die 
Temperature gradients were estimated for the intervals above 
anombelow the upper Devonian im, Fro. |) 2.10.) anc Pftor fehe 
intervals above and below the unconformity surface in Fig. 
2.11. Different gradients were observed for the intervals 
considered. 

Figures 2.12 and 2.13 are temperature/depth plots and 
calculated thermal gradients for the area covering townships 
45, 46 and 47 and ranges 19, 20 and 21. This area falls 
within the disturbed belt, and the plots show the geothermal 
Gracrentse obtained when | the ssstopes sweresacalculatedseior 
Oni cereniumELenologiacal antervals. 

Figures 2.14 and 2.15, tllustrate two rurther such yplots 
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intervals (2) and (4), separated by the upper Devonian 


surface at 3200m. (Townships 51, 52, 53; ranges 19, 20, 21.) 
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Figure Deel) Temperature/depth Lot for, Bathnologaca: 
intervals (2) and (3), separated by the unconformity surface 
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Figure 2.13 Temperature/depth plot for lithological 
intervals (1) and (3), separated by the unconformity surface 


at 3749m. (Townships 45, 46, 47; ranges 19, 20, 21.) 
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Figure Lato Temperature/depth plot for lithological 


intervals (1) and (2) separated by the unconformity 


at 2200m. (Townships 55, 56, 57; ranges 15, 16, ie 
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randes™|.5,- l6 and al7. 

The calculated gradients were used to plot geothermal 
contour maps for the four selected intervals. These gradient 
NepSeaLe Givelsinenics.2 16-26 1/9 em Tnescontour ml atemyalamns 


(227m). 
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Figure 2.16 Geothermal gradient contour map for the 


lithological interval (1): top of the Colorado formation to 


the unconformity surface. 
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Figure LAAT Geothermal gradient contour map for the 
lithological interval (2): top of the Colorado formation to 


the upper Devonian surface. 
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lithological interval (4): the upper Devonian surface to the 


Precambrian. 
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2.1.3 Detailed Temperature Analysis 

Detailed temperature analysis was carried out in order 
to investigate the subsurface temperature distribution. 
Seven levels were considered starting from 2000m to a 
maximum depth of 5000m in steps of 500m. To obtain formation 
Bemperatures Vat weeach | Gepth Vand (sfor = each mote. ene Bx 
township/range areas, the geothermal gradient plots obtained 
BN p62 suhie2ay~WELC— USC .eFOre—each depth, the formation 
temperature which occurred exactly at that depth was chosen. 
Where no actual meaSured temperature value existed at the 
exact depth for a particular 3x3 township/range area, the 
temperature as defined by the linear least squares fit line 
of the temperature/depth plot was taken. 

These temperature values were used to construct contour 
maps of the seven levels considered. The maps are 
Me uSstrateds ie r1lgSmme.. 20MmeO 2026. 

It is informative to determine how the subsurface 
temperatures vary across each of the seven levels. By 
correlating changes im temperature with geological 
formations, deductions can be made regarding the nature of 
the heat flow pattern and its relationship to geology. 

Four profiles were drawi™wacross WEheéw, Seven, Surfaces. 
Three of these profiles are in the SW to NE direction, and 
are labelled AA, BB, and FF. (These profiles are shown in 
Ricemmcec Oe The prorile tabededs pp cutsediagonallyeacnoss 
the entire area while AA and FF are about 50 km from BB 


toward the northwest and the southeast directions 
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Figure 2.20 Contour map of the subsurface temperature 
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Figure 2.21 Contour map of the subsurface temperature 


distribution at 2500m. 
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Figure 2.22 Contour map of the subsurface temperature 


distribution at 3000m. 
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Figure 2.23 Contour map of the subsurface 


Gistribution at 3500m. 
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Figure 2.24 Contour map of the subsurface temperature 


Gistribution at 4000m. 
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Figure 2.25 Contour map of the subsurface temperature 


Gistribution at 4500m. 
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Figure 2.26 Contour map of the subsurface temperature 
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respectively. The limits of the profiles were chosen so that 
only areas where data exist were included. 

The profile labelled PP cuts diagonally across the area 
in the NW-SE direction, and it is perpendicular to the other 
three “profiles. “The elocation of theseaStern Limitmor ithe 
disturbed belt 1svalso vindicated om Figs. 2.20 to 2e@oe 

AYong= cach ot the profiles . Jarid fore eachmmcr the 
surfaces, temperature and formation information were used to 
construct contour plots on vertical planes corresponding to 
each prone. Figures: 2. Ziseec. 20102 3 land 2.33.5 shows “ehnese 
COnCOUMmmMaDS mand. eG. | 2626, 2. o0ME 2 oc and) = 2o dees how 
profiles of the three surfaces mapped in 2.1.2: 

(A) The top of the Colorado formation. 

(B) The surface of the Unconformity. 

(C) The upper Devonian suPtaties 

Since the temperatures used here are not corrected in 
any ner an attempt was made to determine the expected 
error. The Horner-type correction (Fertl & Wichmann, 1977) 
waS applied to temperatures for which there was sufficient 
information. An example of the temperature stabilization, 
and graphical method used to find the equilibrium formation 
Peliperarures,. 15 Shown 1ne@F190.m2.s0 hi GuLe @2. soasgi Ves mabe 
versus time from when circulation ceased to the time at 
which the temperature measurement was made. Figure 2.35b 
gives temperature versus AT/(AT+T) where T is the 
Circulation time and AT is the time from when circulation 


stoppedmeuntil (the. partuwculare measunement eeveasqrmade. The 
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Figure 2.27 Temperature contour map in the vertical plane 


along profile (AA). Temperature in degrees centigrade. 
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Figure 2.28 Topography of the surfaces (A) top of the 


Colorado (B) unconformity surface (C) upper Devonian surface 


along profile (AA). 
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Figure 2.29 Temperature contour map in the vertical plane 


along profile (BB). Temperature in degrees centigrade. 
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Figure 2.30 Topography Of meatne 


Colorado (B) unconformity surface 


along profile (BB). 
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Figure 2.31 Temperature contour map in athe geverticel wa plane 


along profile (FF). Temperature in degrees centigrade. 
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Figure 2.32 Topography of the Surtaces (A) §tOpmOle sue 
Colorado (B) unconformity surface (C) upper Devonian surface 
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Colorado (B) unconformity surface (C) upper Devonian surface 
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Figure 2.35 An example of the Horner-type correction. (From 


Lam et al. 
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equilibrium formation temperature is then obtained by 
extrapolating a least squares fitted line to the vertical 
axis which corresponds to infinite elapsed time. 

This procedure was applied to the wells for which there 
was sufflesent Jinformation, (and the values weresuseduto 
generate the corrected temperature/depth plot which is shown 
Mime ug. 62.36. “The (least “squares fitted! lane sform ene 
corrected data gives a greater slope than that obtained from 
the uncorrected data. To determine the temperature change to 
be expected at any depth, the uncorrected and corrected 
temperature/depth plots were overlaid and the indicated 
correction for the temperature as a function of depth was 
determined. 

Table 2.1 gives the temperature correction that can be 
expected for each depth considered. From this table it can 
bemingerred™that equ: libraumetemperatures, sare shigher = than 
measured. temperatures. Also, the table indicates that the 
temperature differences between corrected and uncorrected 
values increases with depth. In this regard, the calculated 
temperature gradients may be viewed as minimum values. Also, 
from the calculated average relative errors, it is found 


that the uncorrected temperatures are about 12.5% lower than 


the equilibrium formation temperatures. 
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temperature values. (From Jones et al, 1983). 
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Tables? > 1 


Temperature correction as function of depth based on 


Horner temperature correction. 


a eS ee en ee ee 
Depth(m) Corrected Uncorrected +AT°C ¥ Error 
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2000.0 63.0 57.0 6.0 T0r 5 
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2.1.4 Thermal Conductivity Estimates 

In order  Stoesestimate fterrestrial Sheat et lowesiteeis 
necesSary to know the thermal conductivity of the strata 
through which the heat flows as well as the temperature 
gradient across the interval of interest. 

Manyesttata Or difierentimock) —ty pes occ Ursa nEEthe 
sedimentary column in the west central Alberta region. Thus, 
aneetiectivetthermal conductivity for. thegrcolumn | must. ibe 
estimated for each interval across which Heat flow is toebe 
computed. 

BE gcthermalyconductivisy Vvaluesuk,, Ko sakaeieas. pakeecan 
be assigned to discrete layers of.thicknesses L,, Lz, 3, 

-, Ln» respectively, then the effective conductivity for 


the total interval is given by 


Keely / th G/ Ky) | a) 
and 1S equivalent to the weighted mean thermal conductivity 
across the total thickness. 

Jones et al, (1983) measured thermal conductivities of 
936 samples from 48 wells in the study area. The locations 
of the wells are shown in Fig. 2.37. The samples are ‘from 
different formations. The meaSurementS were made on core 
plugs in a divided-bar apparatus. 

In order to obtain effective conductivity estimates, 
these measured values were investigated in detail. 


Initially, the samples were grouped as to rock type, and the 
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Figure 2.37 Locations of wells from which measurements of 


thermal conductivity were made on core plugs. 
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mean thermal conductivities and the standard deviations were 
calculated. The results are shown in the GrapneOtelag. 2.50. 
The number of samples of each rock type is also indicated in 
this figure. Dolomites, limestones, siltstones, sandstones, 
dolomite-limestones, shales, and shaly-siltstones constitute 
about 89% of the rock samples. Among these major rock types, 
the dolomites and sandstones have thermal conductivites a 
little greater than 3.0W/m/K. The limestones are of 
generally lower thermal conductivity. 

Thermal conductivity values representative of each of 
the uri hological intervals, rdéscribed) in Section 2.3 were 
SOUGht.fl0 [éthis* jegard,4 the? thermal™ conductivity as a 
function of depth was investigated for each lithological 
ineer val. sscatcer PlOtSaanershown HnOFigs4s 2539-242. Also, 
these figures show average conductivity values together with 
Standard deviations calculated from individual values for 
all ene oek types in each interval. These diagrams show in 
a general way the depths of the formations from which the 
samples came. Based on these figures, there does not appear 
to be a direct relationship between thermal conductivity and 
depth. 

Pipaddraron,. the relationship between porosity and 
depth was investigated by plotting the measured porosities 
6=t the samples as a function of depth for all samples. | The 


Lesuleingus plot 91S sillustravedmin h1¢ee2.437eand sindicaces 


that there is a decrease of porosity with depth. 
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Figure 2.38 Graph of conductivity 
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conductivity as a function of depth for lithological 


interval (4), the upper Devonian surface to the Precambrian. 
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With these general results in mind, and to better 
understand thes rthermal ‘conductivity distsi bution each 
lithological interval was considered separately and graphs 
of thermal conductivities and calculated standard deviations 
of the measurements for the main rock types in each interval 
were constructed. These graphs are given in Figs. 2.44-2.47 

Comparison of Figs. 2.44 and 2.45 indicates that the 
Siltstones from the Colorado to the unconformity surface 
haveghirgher jhermal conductiuity thanethosesftrom the top) of 
them Coloradom ComithesuppermDevonian This is also true for 
sandstones and shaly-siltstones. For these three rock types, 
there iS an apparent decrease in thermal conductivity with 
depth for the samples from the interval between the 
unconformity surface and the Devonian. 

On the other hand, comparison of Figs. 2.46 and 2.47 
shows that the thermal conductivity of the dolomites is 
slightly greater at greater depths. This is also apparent 
for dolomite-limestones, though the limestone conductivity 
itself does not differ for the two intervals considered in 
these figures. 

Further efforts were made to. find representative 
thermal conductivity values for the lithological intervals 
by grouping the measured values for each interval into 
Various conductivity intervals without regard to rock types. 
All measured values were assigned to thermal Conductivity 
intervals of 0.20W/m/K wide over the whole range of 
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2.46 Graph of thermal conductivity for the main rock 


Figure 


types within the lithological interval (3), the unconformity 


Surface to the Precambrian. 
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Figure 2.47 Graph of thermal conductivity for the main 
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histograms was constructed for the lithological intervals, 
and these histograms are given in Fiqsie 92e0 2 7 iuens 
before, the mean and standard deviation for each interval 
were calculated from the individual values for all the rock 
types in the interval. 

The mean eonductivities calculated by this) last 
approach were taken as the representative thermal 
conductivity values for the lithological intervals. These 
values, together with the main rock types, are summarized in 
Table 2.24 

PaplLe =2 309: Vesmaene “CONdUCLI Vi LY = Vellccmwnie@huEwere 
taken, Skor the  ditterent rock types ‘for “the effective 
conductivity estimates. The sources of data are also 
indicated in this table. The estimated thermal 
conductivities are lower than the average thermal 
conductivities based on measured values. The results, based 
on effective conductivity estimates, take into account the 
proportion of the sedimentary column occupied by each rock 
type, and therefore such estimates are probably better than 
those based on measured values from selected formations 
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Table 2.2 


Representative thermal conductivity fOr four 


lithological intervals based on measured rock samples. 


LITHOLOGICAL GEOLOGICAL 


INTERVAL FORMATIONS 
Colorado to Cretaceous 
Unconformity 
Surface 
Colorado to Cretaceous 
Upper Devonian Jurassic 

Triassic 
Mississippian 
Unconformity Jurassic 
Surface to Triassic 
Precambrian Mississippian 
Devonian 
' Precambrian 


Upper Devonian Devonian 
to Precambrian Precambrian 
& 


MAIN ROCK TYPES MEAN K 
(W/m/K) 


Sandstone 2290 
Shale 

Siltstone 
Shaly-Siltstone 


Sandstone 2 
Siltstone 

Shale 

Shaly-Siltstone 
Limestone 


Limestone Qed, 
Dolomite 

Dolomite- 

Limestone 


Dolomite 2.98 
Limestone 

Dolomite- 

Limestone 


STD 
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2.1.5 Calculation of Conductive Heat Flow 
Conductive heat flow can be calculated from the 


Steady-state heat flow equation 


Q=KaT/dz Zo 
where Q=vertical heat flow, K=thermal Conductivity: 
and dT/dz=temperature gradient. 

To obtain heat flow estimates, average geothermal 
Gradients were calculated forf™each of ™/the lithological 
intervals considered and these were combined with the 
estimated thermal conductivities. The average geothermal 
gradients were calculated for each interval from the 3x3 
township/range values. 

Table 2.5 summarizes the thermal conductivities and 
average gradients for the intervals considered and gives the 
calculated heat flows for these intervals. Also, an error 
estimate is included COG the conductive heat flow 


calculation. This error estimate iS given by 


CO. =0/7 yn 2s 


where 
gn=Standard error of estimate, 
G=Standara deviation, 


n=sample size. 
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Heat flow estimates for the four lithological intervals 


based on the statistical averages of measured conductivities 


and gradT. 


EEE 


LITHOLOGICAL 
INTERVAL 


ERROR 
% 


————————— eee 


Colorado to 
Unconformity 
Surface 


Colorado to 
Upper Devonian 


Unconformity 
Surface to 
Precambrian 


Upper Devonian 
to Precambrian 


MEAN K MEAN HEAT 

(W/m/K) GRADIENT FLOW 

(°C/km) mW /m? 

2.90+1.11 23.86+4.89 69.19 
2.70+1.10 25.70+3.70 69.37 
2774414 30.53+6.64 84.57 
2.95+1.19 31289*8599 95203 
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Table 2.5 indicates a significant difference in heat 
flow above the unconformity Surface compared with that 
below. Also, within the Paleozoic itself, there is a change 
in heat flow with depth as indicated by the results from the 
overlapping intervals from the unconformity surface to the 
Precambrian and the upper Devonian to the Precambrain. 

For comparison, conductive heat flow was calculated 
from estimated effective conductivity values based on 
equation 2.1 and the heat conductivity values from Table 
7 set on 

The calculated heat flow values are given in Table 2.6. 
The results show that there is an increase in heat flow with 
Cepel eDubetOn each Withological s intervaleetne svaluesmare 
considerably lower than those of Table 2.5 due to lower 


average conductivities calculated using equation 2.1. 
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Table 2.6 
Heat flow estimates for the four lithological intervals 
based on calculated effective thermal conductivities and 


average gradT. 


Stratigraphic Average Average Heat 
Interval Conductivity grad T. Flow 
W/m/K °C/Km mWm~ ? 
re eee 
Ground Surface 2.46 = = 


to top of Colorado 


Top of Colorado oe 23,9 55.6 
to top of Paleozoic 


Paleozoic to top 2.39 2527 61.4 
of Devonian 


Top of Paleozoic 2a 30.5 76.6 
to Precambrian 


Top of Devonian 2.49 S469 79.4 
to Precambrian 


3. DISCUSSION AND CONCLUSIONS 


3.1 General Discussion of Results 

The Paleozoic and the Precambrian surface maps of Figs. 
2.4 and 2.5 indicate that the sediments dip in an 
approximately southwest direction toward the eastern limit 
of the disturbed belt. This is further confirmed by the 
prckvles@omh FigsA 2226iie2.30itcand 62.32 sandUimeans  echatmbine 
overlying upper and lower Cretaceous Mesozoic formations 
become thicker toward the southwest. These southwestward 
dipping sedimentary formations extend relatively smoothly 
eastward well beyond the study area. Beyond the disturbed 
belt boundary to the west, toward the Rocky Mountains, the 
sediments are faulted and folded. 

The geothermal gradient map of Fig 2.9 indicates that a 
high geothermal gradient occurs ‘along the edge of the 
disturbed belt, with a regional trend that strikes in the 
northwest to southeast direction. Approximately 
perpendicular to this are three northeasterly trending high 
gradient zones which are indicated by arrows in Fig. 2.9. 

It is apparent that the high gradient geothermal 
anomaly delineated by Lam et al. (1982) and striking in the 
Southwest-northeast direction is part of a larger high 
G@eadrent zone which borders the disturbed belt. 

Thre Mageothe rmall scqradienit fémapsmiicratthe wilrcholed reds 
intervals from the unconformity to the Precambrian surface 
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there are a number of isolated high and low gradient areas. 
simmlarhy, heligshis2 dcadandk2.da7vedo not undrcateua regional 
geothermal gradient distribution that is consistent for the 
two intervals. 

To better understand the thermal regime, the subsurface 
temperature distributions at 500m intervals, as illustrated 
in Figs. 2.20 to 2.26, were considered. These figures show 
that there is a gradual increase in subsurface temperature 
bQOwandreheenortheastuat saowgivent depthhoringnadditionmieut 
pHofilehe(HE)tvofhedFigl; 2¢3treiserconsidered, sathere oistan 
increase of temperature at any given level in the direction 
from the southwest toward the northeast. This implies that a 
horizontal geothermal gradient exists, and the magnitude of 
this gradient is greater near the limit of the disturbed 
bekt. 

ProfilevalAA)neofcFig: t2<27 .showswaesimiulaceresultpebut 
the effect is less pronounced at depths below the top of the 
upper Devonian ene ences Siamalarly, profile) ( BByroterigqat2=29 
shows such a result, but the effect seems to disappear at 
depths below the upper Devonian formations. 

HiaorofhaleulaPPimeofe Big. R2nesn “whachguctnsedininthe 
northwest-southeast direction, 1S considered, there is 
little change in temperature along the profile for any given 
depth. 

By mapping geological surface information onto the maps 
illustrating temperatures for profile (FF), the temperatures 


@kongishe various surfaces can be determined. | From these, 
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horizontal gradients were estimated on the basis of a 200m 
grid. The horizontal gradients along the upper Devonian, the 
unconformity, and Colorado surfaces on profile (FF) are 
STOW Dee) ete. eS one 

EE Of a aii ele eee See OT that the horizontal 
temperature gradient is greater near the edge of the 
disturbed belt, and there is a general decrease in the 
horizontal gradient away from the disturbed belt. 
Purthermore,;”the figure shows that the horizontal) gradient 
near the disturbed belt is greater along the surface of the 
upper Devonian than along the unconformity surface and much 
Yarger than ‘along the Colorado. 

The existence of a horizontal temperature gradient 
indicates that there is lateral heat flow and this may be 
associated with fluid motion. This effect can be due to the 
different hydrodynamic situation between the area close to 
the disturbed belt where most of the downward gravity 
imposed water flow occurs and the shallower part of the 
basin where fluid tends to flow laterally. This implication 
of heat transport by fluid flow is strengthened by results 
from model calculations by Dr. Rahman and included in a 
paper in preparation on results from this area (Jones et al, 
1983). The model considers the influence of gravity induced 
water flow on the temperature as a function of depth. Such a 
Situation, was. first described agrox one dimension by 


Bredelhoeft and Papadopulus (1965) and is represented by the 


equation: 
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4 devonian surface 


X unconformity surface 


© colorado surface 


LATERAL TEMPERATURE GRADIENT °C /KM 


0 20 40 60 SOmeeOOw uiZ0nu 40 
DISTANCE km 


Figure 3.1 The horizontal gradients along the three 


geological surfaces shown Hine spy Zajc! Weyer etetepeMiWe AMS) oh 


ial tele = Panes is 


os wn nag ae 
; ihe : . 


oF 


T(z)=(Ty-To) (1-EXP(-P,z/H)/(1-EXP(-P,))+T, a 


where P, is the Peclet number for vertical flow 


(ie, the ratio of convective flow to conductive flow) 


p. =pCr Shh 0. 


where 


T(z)= temperature at depth z, 

z=depth, 

Ty=temperature at depth z=H, 

To=temperature at depth z=0, 

p =density, 

Cesheat capacity, 

Ac= heat capacity of solid-fluid system and 

Qw=groundwater flow per unit crossectional area. 

The value A./pxC¢ generally equals about 10° ‘m’sec™'. 
Temperature versus depth changes evaluated on the basis of 
equation (3.1) for To=5°C and Ty =150°C and different values 
of downward groundwater flow velocity are given in Fig. 3.2. 
The result shows that even a downward water velocity of 
10-'*m/sec will produce an increase of thermal gradient with 
depth and downward velocities of 0.5x10°°’m/sec and 10° ’m/sec 
give even greater gradient increases With depth. Such 
velocities can explain the temperature gradient increase 
with depth in the area and the discrepancy between the 


average heat flow for the shallower and the deeper 
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Figure 3.2 Temperature dependence with depth in the interval 
between ground level (5 degrees centigrade) and 5km (150 


degrees centigrade) due to downward water movement with 


aifferent velocities. 
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Stratigraphic intervals. 

Tables 2.5 and 2.6 show that there are discontinuities 
in the heat flow at various surfaces. These results indicate 
that lateral fluid movement within the POA S formations 
Of the Mississippian, upper Devonian and the lower 
Cretaceous @plays” an important role in modifying the 
hydrodynamic situation and heat flow pattern: 

A schematic model of the hydrodynamic situation in the 
Scud Vamee (Ca pe USieeG LV.Clie leit G03 oo ae emMOdeimsuggest Sut nat 
Surface water percolates downward in the Rocky Mountains and 
foothills recharge area. The water then moves’ under 
hydraulic pressure along the permeable sediments of the 
lower Cretaceous, MissisSippian and upper Devonian to the 
discharge area to the east and the northeast. 

Such a model appears to satisfy the observations from 
the data, and it does not appear that some other source of 
heat, such aS a magma intrusion or upwelling of partial melt 
occurs. Deep drilling and seismic reflection results (Price 
and Mountjoy, 1970) indicate that the Cambrian basement is 


unbroken. 
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Figure 3.3 Hydrodynamic representation of the study area 


(after Hitchon, 1982). 
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3.2 Conclusions 


The following conclusions can be reached from the work: 
(1) Temperature gradients and heat flow increase with depth 
ine the “study “area. This is especially eoicent in the high 
relief and the high hydraulic head region of the disturbed 
belt and the deeper part of the basin. 
(2) Lateral temperature gradients occur along the main 
stratigraphic surfaces. 
(3) The increase in heat flow with depth and the lateral 
temperature gradients can be explained by the hydrodynamic 
Situation which is characterized by downward gravity imposed 
water movement to great depths and probably upward and 
lateral water movement in the deep permeable strata and 
Ehrough faults. 
(4) Petroleum bottom-hole temperature data analyses in 
sedimentary basins are useful for detecting major heat flow 
features related to the basin hydrodynamics and are a 
necessary first step to determine the significance of 
detailed individual heat flow meaSurements. 
roy Lt LS Snot Pesutficient Peto. idraw Tconciusronsesabout 
terrestrial heat flow from a small number of even extremely 
accurate temperature meaSurements which come from depth 
MLE hValGmewhich «i NClude mon yam Oat ta Ones sune sedimentary 
section. The heat flow pattern within the strata can be 
significantly disturbed by the hydrodynamic situation. 
(6) Steady-state “background” heat flow in the region can 


only be determined from careful measurements below the 
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sedimentary strata, ie, in the deeply buried Precambrian 
basement formations which may be less influenced by the 


water flow system. 


iadansexd pa 


7 


dooasettabs 


pene = ith 


BIBLIOGRAPHY 


American Association of Petroleum Geologists, PSMA SA 
Bottom-hole temperature data file. University of 
Oklahoma. 

Ang lamyisgrQiiM’A@gande,.Beck-dOa. =. 1965. Regional Heat Flow 
Pattennein® Western (Canada. Can. “U. “Farth sSoie, a2: 
WLO-RO2Z. 

Ballang wn. ,VWktret@ansen, (I86hingpoulseny D. K., and Saxon, 
S. 1980. Proceedings of the II International Seminar on 


the Results of EC Geothermal Energy Research. Strasburg, 
March, 1980, Reidel Publishing Company. 


Puroip FMR. Roy, FR. RepandgDeckeGgmerheRmal 96S meteanrer] ow 
and Thermal History in New England and New York, In 
Studies of Appalachian Geology, ed. E-en Zen, pp. 
437-451 Interscience, New York. 


Bodvarson, G. 1969. On the Temperature of Water Flowing 
Through Fractures, Ue “SOFA Geophys= i kes, “if pp. 
96799925 


BOUCemeMAIVaH. JeP Raw, 1 982F%enihneotinterlorn ofesahke Earth its 
SEPUCEUrE. *COASLITULIONW and “EvoluGion» SEdward sAnncid 
(Publishers) Limited, London. pp. 266-295. 


Bredehoeft,) J. D., and Papadopulus, T. S. 1965. Rates of 
Vert ical Ground Water Movement and Estimates from 
Earth’s Thermal Profile. Water Resources Research, 4: 


S273 oF 

BRGM- 1976 Report. Potent ial Geothermique du _ Bassin 
Parissien, Bureau de Recherches Geologique et Minieres, 
Paris. 


pumlardweC.. Be 939. Hear Fiow li SOUGTIEALL ICa sy PrOC mm Roy. 
Soc. London, A173(995) pp. 474-502. 


Burwash, R. A., 1979. Uranium and Thorium in the Precambrian 
Basement of Western Canada. II. Petrologic and Tectonic 
Gonnrorse Canttun-Earth Sempre? 22464: 


Bouwash,eR. Av, cand Cumming, CG... 1976. 8UranlUmeaid ey) non in 
in the Precambrian Basement of Western Canada. Can. J. 
BanthaScHetererl234-293t 


Canva lion M. De, sPUrWOKO) 85s, Thamrin, M.;.and Vacquier, vie 
1980. Terrestrial Heat Flow In The Tertiary Basins OF 
Central Sumatra, Tectonophysics,69: pp. 163-188: 


Elder, J.W. 1981. Thermal Systems. Academic Press, Bins Ax 


oy. 


Ne oe 
in -ythenew - | 
¥ - ; | 

- alas bre Sons ae rear 

‘aie at SD Mined 2a “i ee seo. 1. : ikea ii . 

| vanish = tian ;s : 

4 Sh yee 2 ..A9@ ns ‘ i. 7 

he lage. | esd a bens : a egned inlay — 

_ 

= 

7 


iene Ronee Gp va Veune eet = 


viol" te 2eet A 

ol ,ange wed Yas hit pm aa a 

G9 .7@S na-a be haa 
: ost saa Sazeent Teb-teER : 


. = 
siuol’. Seta. to sar ‘oad © .293987 2 .roaavbol i 
; - BaP PNAC A : 
; 
_ 
7 


ie seer -r8er wan 
etl Area oft Fe Bak: | <9 ae M a ee a 
aah S otal. er ees _—- : 

| mm oh ad el aint (erat nak an 
40 e@Ysh .dg2! .2 .P utoqetaget Bie. .d .o .2ipotsbere a 
aN 2eaTeMtes. ‘bes Vee ey Pr on +e 
:b ,dotesget eetuyews® setew . Si) Vesnnver th Ye ee 


Meesa. A 
eansinit: 1S 


98 


Dome 1-36, 


Evans, T. R., and Coleman, N. C. 1979. North Sea Geothermal 
Gradients, Nature, 247: pp. 28-30. 


Fertl, W.M., and Wichmann, P. A. 1977. How to Determine 


Static Bottom-hole Temperatures from Well Log Data. 
World Oil, 184: pp. 105-106. 


Gablew Ra! 9/9. me part Or a Geothennaleni tux Map of France 
pp.179-186.°In Terrestrial Heat Flow OF Europe Feck, 
Rybach, V Cermak eds.) Springer, 1979. 


Garland, G. D.,and Lennox, D. H. 1962. Heat Flow in Western 
Canada. Geophys. J., 6: pp. 245-262. 


Harper, M. L., 1971. Approximate Geothermal Gradients in 
North Sea Basins, Nature 230: 235-236. 


Hutcheon 8.) 1982. 8hiulid flow in Western = Canada Sedimentary 
Basin. A, Effect of Topography. Water Resources 
Researchy5: 186-795" 


inghantpesMiey Re, “SBinghanyD. Kee and iGough) DBD. el, S83. 
Magnetovariational Study of a Geothermal Anomaly. 
Geophysh. JG Roy avAstroenmteSoc #u72. (Uintpress) 


Jessop, FApeMsmetobDar re GMatA., Pand’Selatenmeu. Go 1976.) ihe 
World Heat Fiow Data Collection. Can. Earth Physics 
BoanGhigeoy 612 OBO. 


Jessop, A. 1976. Geothermal Energy from Sedimentary BaSins. 
Can. Earth Physics. Geothermal Series 8; 10pp; Ottawa. 


JONES me fa WwW. tar Rahnian pp MayeMajorowic?, fil. TAS, Sand Lam, H. L. 
BOGS. PEREPORT® “ABGSEUdG) Ser Thermal Properties of 
Sedimentary Rocks of Western Canada" (UP-A-269-PHASE II 
) 


Jones} FORW., Kushtgqbor , MCCAyLamy RH leu0e Magiorowiczecy. Ae 
and Rahman, M. 1983. Estimates of Terrestrial Thermal 
Gradients and Heat Flow Variations with Depth in the 
Hinton-Edson Area of the Alberta Basin Derived from 
Petroleum Botton-Hole Temperature Data. Submitted to 
Geophys. Prosp. 


Kappelmeyery (60. Nand Haenel /% R#OrTST4N) UGeothermicseewi tn 
Special References to Application. Gebruder Borntraegen, 
BeclimsestubidarO, poerader envy 2 


Dachenbrtich, wAvl Hesand Sass. HIBSLITE. Heat Criowr ine the 
United States and the Thermal Regime of the Crust. In: 
FRérEarth’ scCrostitud® 'Grptéaccekt Ed. |AmereGeopnys. 
Monograph 20, No.8, pp. 626-675. AGU, Washinton D.C. 


we es 


pein gle wok, ay BF al nt 7 oe. ame 7 te a al 
ATG ¥ sfori-mottes eT) 


eore= Ao ge xu t 
7 e econud 


yn R VSR { wad Roane. 
zanrTuoaen. TS78W 


VV enter Tah 


3atq ff 


ayy 3, ener . BD) VG staat thy ote 
anieyad Monel | se ees 


io pivart i ornswse 


.2- SR ,med bre Mati cesane ry 
10 Bet yneeaw | 
i sees cna 


a - 
cM. 
seat ‘bl 


Al segivere get 
Venrarit \eletas 


7 aad a re 2a - a a i a 


99 


baliponeul. , Jones, F. W. and Lambert, C. A. 1982. Geothermal 
Gradients in the Hinton area of West central Alberta. 
Can. Sa7tombRarth=Ser. 192/755-7766. 


Meyorewlez, “OMEASY, Jones,.F #lawieeland Lam, @H@ DBererossp., 
Terrestrial Heat Fiow in Relation to Hydrodynamic in 
Alberta Basin Presented at the Proc. of IUGG Hamburg 
Conference. 


Ma jOPowlcz rd Ay! ConesHintorwes, andeyLamee He. L.  1983a. 
Comment on a Magnetovar iat ional Study of a Geothermal 
Anomaly by M. R. Ingham, D. K. Bingham and D .I. Gough. 
pears LorUcpub lacaniononi natheaGeophvs. aoe RABASErO! 

OGS, ; 


MayOrowloz sd < A., and Jessop, A. M. 1981a. Regional Heat 
Flow patterns in the Western Canadian Sedimentary Basin. 
Tectonophysics 74: 209-238. 


MajoOmowlcz | tigeA., and Jessop, A. M. 1981b. Present (Heat 
Flow and a Preliminary Paleogeothermal History of 
Central Prairies Basin of Canada, Geothermics, Vol.10, 
Noeenapptesi-938e4 London: 


MeCrossan, RR. G., and Glaister, RR. P. 1970. Eds. Geologica) 
History of Western Canada. Published by Alberta Society 
of Petroleum Geologists, Alberta, Canada. 


O' Brien, PB. J., 1970.Aquifer Transmissivity Distribution as 
Reflected by Overlying Soil Temperature Patterns: Ph.D 
thesis, The Pennsylvania State University. 


Oxburgh, E.R., and Andrews-Speed, C. P. 1981. Temperature, 
Thermal Gradients and Heat Flow in the Southwest North 
Sea. Petroleum Geology of the Continental Shelf of 
North-West Europe. pp.141-151. 


Puiceminke A, and Mountjoy, E.. W. 119702 sGeologicals structure 
OF the Canadian Rocky Mountains between Bow and 
Athabasca Rivers-A Progress Report. In Structure of the 
Southern Canadian Cordillera. Ed. J. 0. Wheeler. The 
Geological Association of Canada, Special Paper) 6, pp. 


Vag Oe 

Reuter er OOLge, sCse wand. (Oval jes hen glo 2. EStimatess of 
Terrestrial Heat Flow in Northern Chihuahua, Mexico, 
Based upon Petroleum Bottom-hole Temperatures. 


Geological Soc. of America Bull. 93: pp. 613-624. 


Reiter, M., and Jessop, A. M. 1983. Estimates Of Terrestrial 
Heat Flow in Offshore Eastern Canada derived from 
Petroleum Bottom-hole Temperature Data. (not published). 


SET vik : 
Vasari iS A Yeon 


Sty & NO. 7 
iouod 57. @ 3G ® ae | pl areinol nf 
poral wy dey notte 78 an 


taal) ipnaloah .st Gtr vary a 
ieee wie. tem Se elect »:: 
me ar ry A Grint Bap y.A 6b isaiwoxoteM ; 
tc «= yoetbe a, igang panes ie 8 Ons 7 


Ot ON, telategethoae 40 Wizad col nis taained a 
M 4 ice saat R= TS 8S a ) 
wigs nee. ,BBe OFeh Ff /deebetolbos oO a taeeovson ae 
else ayueee ‘yh bane ie aberad qrsiesk to yrowald | 
eta? . (SENeeTS esqnetaas ue lasses aouay 7 
as, | neh hols - Ng soviitn’ 


panes , tes9 ind 38 i evra 
Bas ees 1 SAteNs ya betes tea a 
| | a I . ae | he = : : 


sae 


2, &) Pa 


aie nant aS — . 

now sanMeuce = nl i 2 aaneltend 

16 Lait Deprani med ens 30 
gee 


- ye tg ae ; 
-! — Senin 
Prt bert 


100 


Roberts, W. H. 1979. Some Uses of Temperature Data in 
Petroleum Exploration: Paper Presented to Sympossium II, 
Unconventional Methods in Exploration for Petroleum and Gas, 
Dallas, September i3-i4, 1979) 


Rov oh oe DECkCipmbeuRk st, sb lackwe ll G |) Dee mandaepar chr. 
1968. Heat Flow in the United States, JU. of Geophys. 
RES) 73 520/725271,, 


Rybach, L., and Muffer, L. J. 1981. Ed. Geothermal Systems: 
Principles and Case Histories, Wiley, 359 pp. 


Selater as G.) me Jaulpart (Cand Galcon- Deu o6 0m iceneat 
Flow Through Oceanic and Continental Crust and the Heat 
Loss of the Earth. Rev. Geophys. Space Phys., 18: 
269-5) 11k 


Scleter, J. G., and (Parsons, B. 1981. Oceans and Continents: 
Similarities and Differences in the Mechanisms of Heat 
HOSS: aie of GE Oplivy See RES VO ly, CO; mmINO, mb ean DO: 
i eee Paes Hl 


Sie ance Chapman, «Dees aml0ce. wOnu lnernale st Fecuslor 
Groundwater Flow on Regional Systems, submitted to J. Of 
Geophys. Res., May 1982. 


Van Ostrand, C. E., 1939. Temperature Gradients In Problems 
of Petroleum Geology, pp 989-1021, Amer. Assoc. Petrol. 
Geol., Tulsa, Okla, 1939. 


Moles, TEA ila n Eilelell akigetelitwis 8, 8h WER SiiekiMelaiia= elie) 
Thermal History of Queen Charlotte Basin, Can. J. Earth 
Sere 2 OE wales oa1 OS. 


ee . - A 


we 


‘ = 
, See 
pee or a ke % 
: ey : ~ Py 
‘ Dei x's thic* és 7 § ¥ ’ 
; 1s : a ee ee ee 


. 


a 
a 


oO fF “ ! - : 
ri. eSROTC Ree a 


7 


Yo 2158783. lenatyond . 
Jo Bert 


j R 


anstena® oT etnetossa: at 
lente" .2oeeR seg’. TSOn 


dared tage griteetral 


~ “ee = 4 
- $3 | 7 3 
“C2 4 wh : 
bs 7 
—— + = 5 a : 
pare lgne: = UE - 2 — ead Se eRe 
: “<) . a — ; a , 
ea Fs 
5 ~~ = 
7 
— 
i 
- 7 
~ 
Bs 
<= 
sai 
“ agit 
i 
- 2 
= " = : 
4 e 
rh as 
See 7 
= —= =| 
y af na + 
we a 
-_ ~ a = | 4 
hg mined 7 4 — a 
z aie ve ~ 
m ~ — ~ “ = 
. 7 ma ese . 
: i oe : 
—_—. > 7 a) : 
_ a i 
bal —— > - a 7 - 
_ i _ 5 
ee oe, 7 as: 
me =. ei — - a ; oa q 
po Pe as A a 5 a ~ 
aA aan Ce as See ¢ a — . 
4 Sei = = . 
ee ee, eee 
=. i = Ss Le, —_ ae = 
a ea ee Ae : het 


rn Ae ma 
- 


gene yen Sees. As i an 


Table 
te asthe te Hr 
Kistalie aie oe 

eres 


ie 


Hine 
eat Ge i 
nee 


sate 
a 


4 
eWay 
Bare) 


peridiy et 
Heine Mr 
| Tusa 


aU bid li 
ANN 


Eh 
ie 
ibe 
ty 
4) 


aii 


By ale 


Vwhyna haga 
ve 


sinintls Gayl s 
HONE Faactyd pa 
Weal 

he 


ALE 
MR Way 


Hy 
i 
Wig hi 


Pvt ay Yi Bia) 
OA CO 


AT Rts 


a 


ROU 


eT 
a 


